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Abstract We examine the regenerative cutting pro- 
cess by using a single degree of freedom non-smooth 
model with a friction component and a time delay term. 
Instead of the standard Lyapunov exponent calcula- 
tions, we propose a statistical 0-1 test analysis for chaos 
detection. This approach reveals the nature of the cut- 
ting process signaling regular or chaotic dynamics. For 
the investigated deterministic model we are able to show 
a transition from chaotic to regular motion with in- 
creasing cutting speed. For two values of time delay 
showing the different response the results have been 
confirmed by the means of the spectral density and the 
multiscaled entropy. 
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1 Introduction 

A cutting process is a basic machining technology to 
obtain the surface of the assumed parameters. In cer- 
tain working conditions it can be disturbed by chatter 
appearing as unexpected waves on the machined surface 
of a workpicce. The appearance of chatter was noticed 
and described by Taylor in the beginning of 20th cen- 
tury pQ . But the first approaches towards explanations 
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of this phenomenon came about 50 years later through 
the analysis of self-sustained vibrations [2], regenera- 
tive effects [3J, structural dynamics [3J[S], and, finally, 
the dry friction phenomenon [HIE]- Consequently, elim- 
ination and stabilization of the associated oscillations 
have become of high interest in science and technology 
[51IMTU]. The plausible adaptive control concept, based 
on relatively short time series [TT] , has been studied to 
gain deeper understanding. 

Recently, apart from the widely developed chatter 
vibrations chaotic oscillations caused by various system 
nonlinearities were predicted and detected [T2JH3UH1 
UMIMIZUIMiniUSI - The recent technological demand is 
to improve the final surface properties of the workpiece 
and to minimize the production time with higher cut- 
ting speeds [3T]. Thus a better understanding of the 
physical phenomena associated with a cutting process 
becomes necessary |22| . In this paper we will continue 
the work on chaotic instabilities in cutting processes 
proposing the 0-1 test [23 l [24 l [25 l l26 | [27] as a tool identi- 
fying a possible chaotic solution [2"5] . 

This paper is organized as follows. After the present 
introduction (Sec. 1) we describe the model in Sec. 2. 
In Sec. 3 we provide the results of the simulations and 
corresponding power spectral densities (PSD) while in 
Sec. 4 the 0-1 test is applied and, subsequently, the find- 
ings are confirmed by means of the multiscalc entropy 
(Sec. 5). The paper ends with conclusions (Sec. 6). 



2 The model 

A regenerative cutting process may exhibit a wide range 
of complex behavior due to frictional effects [5UT2"] , struc- 
tural nonlinearities and delay dynamics [T71H51I3TJ1 
|3"TU3"2"] . Moreover it may also involve loss of contact be- 
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Table 1 Parameters used in the model 



Fig. 1 Physical model of a regenerative cutting process 1171 . 



Parameter 



Value 



initial cutting depth ho 

frequency of free vibration ojo 

damping coefficient c 

effective mass of the workpiece m 

friction coefficient ci 

chip width w 



lCT :i m 
816 rad/s 
86 Ns/m 
17.2 kg 

1.25 x 10 9 N/m 2 
3.0 x 10 -3 m 



tween the tool and the workpiece. The following equa- 
tions model the regenerative cutting process and the 
mentioned properties. 

After the first pass of the tool, the cutting depth 
can be expressed as 

h(t) = h -y(t) + y(t-T), (1) 

where y(t—r) corresponds to the position of the work- 
piece during the previous pass, and r is the time de- 
lay scaled by the period of revolution of the workpiece 
27r/i?o (Fig. 1). The motion of the workpiece can be 
determined from the model proposed by Stcpan [30j 

ij + 2-fy + (j%y = —sga(v - y)(F y (h) - F y (h )), 
F y (h) = 0(h) Cl wh 3 ^, (2) 

where ljo = ^Jk/m is the frequency of free vibration, vq 
is the feed velocity, and 2j = c/m is the damping coef- 
ficient. F y (h) is the thrust force, which is the horizontal 
component of the cutting force, and m is the effective 
mass of the workpiece. The thrust force F y is based on 
dry friction between the tool and the chip. It is assumed 
to have a power law dependence on the actual cutting 
depth h and to be proportional to the chip width w and 
a friction coefficient c\. &(■) denotes the Heaviside step 
function. The restitution parameter (3 = 0.75 is associ- 
ated with the impact after contact loss, while t~ and t + 
denotes the time instants before and after the impact. 
Substituting Eq. ([T]) into Eqs. we derive a delay dif- 
ferential equation (DDE) for the workpiece motion y(t). 
Plugging its solution into Eq. ([TJ results in the history 
of cutting depth h(t). 



3 Simulation results 

The non-smooth model equations are solved by a sim- 
ple Euler integration scheme. The used parameters [17l 
152] are presented in TablcQ] Furthermore, the feed ve- 
locity vo has been assumed to be fairly large so that 
v >y. Note that, in this case, the system nonlinearities 



are limited to the exponential dependence of the cut- 
ting force on the chip thickness and to the contact loss 
between the tool and workpiece. 

The corresponding time series for two choices of the 
time delay parameter r = 1.8 and 2.1ms are presented 
in Fig. [5] These series have been plotted with points. 
On the first sight one can notice that both solutions 
are complex but the Fig. [5^, shows points grouped in se- 
lected lines while the distribution of time history points 
of Fig. [2}) looks more random. In Fig.[2|3 h reaches nega- 
tive values that signal that the contact between the tool 
and the workpiece is lost. The power spectral densities 
(PSD) of cutting depth S{ui) = 2ir/T \F{h(t)}\ 2 for the 
two chosen delay times (r = 1.8 ms and r = 2.1 ms j3 in- 
dicate a transition from regular to chaotic motion. The 
sharp peaks in Fig. [3] belong to a high-periodic orbit 



J^i-} denotes the Fourier transform. 
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Fig. 2 Time series of cutting depth h(t) for (a) time delay 
r = 1.8ms and (b) r = 2.1 ms (sampling time At = 1 ms with an 
integration step Zli/10 3 =l/xs) indicating regular and chaotic mo- 
tion, respectively. The time series are plotted by sampling points. 
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t = 1 .8 ms 
% = 2.1 ms 



Fig. 3 Power spectral density S(uj) for two chosen delay times. 
A broad band spectral density indicates chaotic / stochastic dy- 
namics whereas sharp peaks imply regular motion. 



growth rate K in the limit of large times 



A' 



M(n) = lim — V [{pj +n -pjf + {q j+n -qj) 2 } , (4) 

iv— >oo 1\ £ — * 

3=1 

ln(M(n) + 1) 



K= lim 



Inn 



(5) 



For almost all values of the constant cq the parameter 
K is approaching asymptotically or 1 for regular or 
chaotic motion, respectively. 

Note, practically, one has to truncate the sums in 
Eq. (gj). Thus we derived K w 0.21 for r = 1.8 ms and 
K « 1.09 for r = 2.1 ms, which supports the first im- 
pression gained from the time series themselves, Fig. [2]i 
and b. Note further that for delay time r = 1.8 ms, 
K decays with increasing n on much smaller values, 



(regular motion) whereas the broad spectrum indicates 
chaotic dynamics. 

Both power spectra are dominated by a main peak. 
In case of regular motion its position belongs to the de- 
lay time t = 1.8 ms while in case of chaotic dynamics 
the time scale belonging to the peak (t p « 2.0 ms) is 
smaller than the delay time r = 2.1ms. This smaller 
value could be a consequence of a tool-workpicce con- 
tact loss. Based on that we take a closer look on other 
measures to characterize the model's dynamics and use 
a 0-1 test for chaos to display a possible transition from 
regular to chaotic motion with increasing delay time r. 



x10 b 



(a) 




x10 6 



4 Application of 0-1 test 

Based on the time series {hj} which is a discretization 
of the solution h(t) of the DDE normalized by its stan- 
dard deviation, we define dimensionless displacements 
in the (p, q)-plane in the following way 23,24 ,28 

n n 

Pn = ^2 hj cos(jco), q n = ^2 hj sin(jco), (3) 

3=0 j=Q 

where Co is a constant. In this way regular dynamics is 
related to a bounded motion while any chaotic dynam- 
ics leads to an unbounded motion in the (p, q)-j>\aric 
[23], see Fig.gk, 

To obtain a quantitative description of the exam- 
ined system we perform calculations of the asymptotic 
properties defined by the total mean square displace- 
ment (MSD) M(n), Fig. go, and finally we obtain the 
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Fig. 4 (a) For regular motion (r = 1.8 ms) trajectories in (p,q) 
stay bounded around the initial point '1'. For chaotic time series 
(r = 2.1 ms) trajectories in p and q coordinates show Brownian 
motion-like behavior, (b) Thus the mean square displacement 
M(n) increases with time. For estimation of p n and q„, Eqs. (0, 
we used Co = 0.7 and for M(n) and K, Eqs. and J5J, the upper 
limits of N, n are N max = 40000, n max = 280000, respectively. 
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Fig. 4 (c) For regular dynamics (r = 1.8 ms) the if -value con- 
verges towards zero with increasing time. For chaotic dynamics 
(t = 2.1 ms) K stays close to 1. 



155] . which improves the convergence of the test (Fig. 
4c) without the consideration of longer time series, to 
find the time delay r leading to chaos (see Fig. [5]). We 
defined a modified square displacement D(n) which ex- 
hibits the same asymptotic growth 



D(n,c ) = M(n,c ) - V osc (n,c ), 



(6) 



where the oscillatory term V osc (n,co) can be expressed 

by 



t t , s 7^12 1 - cos(nco) 

V OS c{n,c ) = E[h\ -— , 

1 - cos(c ) 



(7) 



and E[h] denotes the average of examined time scries 
hi 



0.75 - . 
* 0.5 - 
0.25 - 

* . » 

1.7 1.8 1.9 2 2.1 2.2 2.3 

t [ms] 

Fig. 5 A"-valucs for different delay times r indicate a transition 
from regular to chaotic dynamics in the region of 1.9 ms. 

Fig. 2b, which corroborates the result pointing towards 
regular motion. 

Note that, the parameter cq acts like a frequency in 
a spectral calculation, cp. Eqs. ([3]). If it is badly chosen, 
Co/ At resonates with one frequency of the process dy- 
namics h(t). Such a frequency belongs to a a peak in the 
PSD, Fig. [3] In the 0-1 test regular motion would yield a 
ballistic behavior in the (p, q)-plane and the correspond- 
ing quadratic growth of MSD results in an asymptotic 
growth rate. The disadvantage of the test, its strong 
dependence on the chosen parameter Co, could be over- 
come by a proposed modification. Gottwald and Mel- 
bourne [241126] suggest to take several randomly chosen 
values of Co and compute the median of the belonging 
K- values. Particularly, in Ref. [55] the problems of av- 
eraging over Co as well as sampling the data points are 
discussed extensively. We followed this approach [35] 



1 N max 

E[h] = —— £ Jh, (8) 

max - . 
i=l 

where N max is the number of hi elements. Consequently, 
the oscillatory behavior is subtracted from the MSD 
M(n, Co) and the regression analysis of the linear growth 
of D(n, Co) (Eq. 6) with increasing n is performed using 
the linear correlation coefficient which determines the 
value of K CQ . 

K C0 = -^EM , (9) 
v /var(X)var(D(c )) 

where vectors X=[l, 2, n max ], and D(cq)= [-0(1, Co), 
£>(2,c ), D(n max ,c )}. 

The covariance cov(x, y) and variance var(x), for 
arbitrary vectors x and y of n max elements, and the 
corresponding averages E[x\ and E[y] respectively, are 
defined 

cov(x,y) = — !— ]T (x(n) - E[x])(y(n) - E[y}) 7 

n— 1 

var(x) = cov(x, x). (10) 

Finally, the median is taken of A" Co -values (Eq. 9) 
corresponding to 100 different values of Co € (0, tt). The 
results of K for different delay times r, Fig.JSJ in the 
window between 1.75 ms and 2.3 ms indicates a transi- 
tion from regular to chaotic dynamics with increasing 
delay time in the region of 1.9 ms. 

As a consequence we conclude that in the investi- 
gated window increasing cutting speed leads to a tran- 
sition from chaotic chatter dynamics to regular motion 
with improved surface quality. 
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5 Multiscale entropy 

To characterize the solutions of the DDE, Eqs. ((1) and 
P]). with regard to information production rate and 
complexity, we aim to calculate multiscale entropy (MSE) 
[33]. This method was successfully applied to analyze 
the complexity of biological signals [MII35) . It is suit- 
able for short and noisy time series. As a consequence 
the chosen procedure would be applicable to experi- 
mental data as well. We use an algorithm provided by 
PhysioNet [36]. First we compute coarse-grained time 
series {i^'} using non-overlapping intervals contain- 
ing TV equidistant data points hi, 

4"' 4 £ (") 

t=(j-l)iV+l 

In the next step we calculate sample entropy [37] 
for these coarse-grained time series. Sample entropy is 
the negative of the logarithm of the conditional proba- 
bility that sequences of to consecutive data points x^' = 
(x^ N ', . . . jX^^-i) an d xj. close to each other will 
also be close to each other when one more point is added 
to them. Hence it is estimated as follows 

^W)=-ln%±A (12) 
Um ir) 

where Um\r) represents the relative frequency that a 
vector is close to a vector X^ (i ^ j). Close to 
each other in the sense that their infinity norm distance 
is less than e = ra. By a we denote the standard de- 
viation of the data. In the limit of m — > oo and r — > 
sample entropy is equivalent to order-2 Renyi entropy 
Ki and is suitable to characterize the system's dynam- 
ics [35]. For independent variables {£} the entropy fol- 
lows from S { /\m,r) = - In P(\^ N) - \ < e) and 
is independent of word length to. For Gaussian white 
noise (GWN) the coarse-grained time series is known 
to be Gaussian distributed too. For small e this yields 
S^\m,r) w -^^/(VttctW)]. Using that the stan- 
dard deviation of the coarse-grained time series erW 
decreases with 1 / s/N leads to following expression 

s£°(m,r) « -ln(rW^) (r -> 0). (13) 

To clear up the characteristics of the cutting pro- 
cess, we look at MSE depending on box size r for the 
two chosen delay times, Fig. [6] For regular motion we 
expect the entropy to approach zero with decreasing 
r. This is observed for the time series with delay time 
r = 1.8 ms. For chaotic dynamics the entropy should 
stay finite, observed for r = 2.1ms. For the sake of 
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Fig. 6 Multiscale entropy Se depending on scale factor N and 
box size r. For regular motion we expect the entropy to approach 
zero with decreasing r. For chaotic dynamics the entropy should 
stay finite. Since Se is not decreasing with scale factor signifi- 
cantly, it seems there is no characteristic time scale present in the 
data. 

completeness, it should be mentioned that in the case 
of stochastic dynamics the entropy would diverge with 
decreasing spatial resolution r, cp. Eg. (|13|) . In Fig. [6] 
and [7] we further analyze the scale factor dependence 
of MSE. The entropic measure is always larger for the 
chaotic time series since it is the more complex one. 
MSE for small scale factor, Fig.[7Jt, indicates that there 
is no characteristic time scale, comparable to l//-noise 
[33] ■ But for larger scale factors MSE is decaying com- 
parable to Gaussian white noise, Fig.[7}}. Thus, even in 
the chaotic case, there exists a characteristic time scale 
which is close to the delay time. The frequencies domi- 
nating S(u>), Fig. [31 are also present in S'^ V ^(to, r). They 
belong to minima in Fig. [7] We learn coarse-graining of 
the data over multiple of time scale t p belonging to 
structures in the cutting process dynamics leads to less 
complex time series and contains less information. 



6 Conclusions and last remarks 

Concluding, the 0-1 test differentiates between the two 
types of motion. Depending on the chosen delay time 
for the investigated DDE, Eqs. (UJ and (0), regular or 
chaotic motion is observed and a transition from chaotic 
to regular motion is detected with increasing cutting 
speed. The nature of solutions has been also confirmed 
by the corresponding power spectral densities and mul- 
tiscale entropies. The latter reveals more insights into 
the process dynamics but is of much higher computa- 
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5 10 15 20 25 
(b) NAt [ms] 

Fig. 7 Multiscalc entropy Se for fixed m = 2 and r = 0.05 
depending on scale factor N. The time series with delay time r = 
2.1 ms seems to be more complex than the time series with delay 
time t = 1.8 ms since its entropy is larger, a) To gain a higher 
scale factor resolution we recorded also time series with smaller 
sampling time At = 0.1ms. b) The existence of a characteristic 
time scale which is in order of magnitude of the delay time is 
indicated by the decay of entropy with increasing scale factor N. 
The gray line represents GWN with r = 0.25. 

tional cost than the 0-1 test and the spectral calcula- 
tions. 

The 0-1 test appeared to be relatively simple and, 
consequently, useful for systems with delay and discon- 
tinuities. A huge advantage of the test is its low compu- 
tational effort and the possibility to compute it "on the 
fly" while the data is still growing. One of the useful 
aspects of the 0-1 test is that the result can be plotted 
against the parameter r. 

The presented method gives a quantitative criterion 
for chaos similar to the maximum Lyapunov exponent. 

As demonstrated by Falconer el al. ^?5\ and Krese 
and Govekar [33] the method can be used on experimen- 



tal data as well. Unfortunately in case of the cutting 
process experimental data are often characterized by a 
relatively high level of noise [20]. In the examined sys- 
tem, we waived the possibility of additive noise. It was 
shown that the 0-1 test could be applied on dynamical 
systems with additive noise and a good signal to noise 
ratio [23] ■ 
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